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SUMMARY 


Aerodynamic  research  relating  to  modern  helicopters  includes  the  study  of  three- 
dimensional,  unsteady,  nonlinear  flow  fields.  A  selective  review  is  made  of  some  of 
the  phenomenon  that  hamper  the  development  of  satisfactory  engineering  prediction 
techniques,  but  which  provides  a  rich  source  of  research  opportunities:  flow  separa¬ 
tions,  compressibility  effects,  complex  vortical  wakes,  and  aerodynamic  interference 
between  components.  Several  examples  of  work  in  progress  are  given,  including 
dynamic  stall  alleviation,  the  development  of  computational  methods  for  transonic 
flow,  rotor-wake  predictions,  and  blade-vortex  interactions. 


I.  INTRODUCTION 


The  flow  fields  of  helicopters  provide  some  of  the  most  complex  challenges  to  be 
found  in  the  field  of  applied  aerodynamics.  For  moderate  flight  conditions,  a  combi¬ 
nation  of  practical  experience,  empiricism,  static  airfoil  characteristics,  and  lin¬ 
ear  theory  are  usually  adequate  to  estimate  the  overall  performance  and  to  assess  the 
relative  merits  of  various  configuration  changes.  However,  it  is  the  boundaries  of 
the  flight  envelope,  which  are  typically  set  by  vibrations,  excessive  power  require¬ 
ments,  aeroelastic  instabilities,  and/or  adverse  handling  characteristics,  that  often 
determine  the  operational  success  of  modern  helicopters.  Expanding  these  boundaries, 
in  turn,  generally  involves  improving  aerodynamic  characteristics  that  are  limited  by 
very  complex  and  nonlinear  phenomena,  such  as  boundary  layer  separation,  and  aerody¬ 
namic  interference  between  components  of  the  helicopter.  These  sundry  phenomena  are 
indicated  schematically  in  figure  1. 

This  paper  describes  a  few  of  the  special  research  problems  and  opportunities 
that  will  be  important  if  major  improvements  in  helicopter  characteristics  are  to  be 
made  in  the  future.  The  list  of  topics  described  in  the  following  sections  and  in 
the  references  are  only  intended  to  be  representative,  and  perhaps  provocative, 
rather  than  complete  or  exhaustive.  Furthermore,  the  material  presented  tends  to 
reflect  the  author's  personal  interests,  biases,  background,  and  research  activities. 
In  each  of  the  three  general  categories  of  viscous  flows,  transonic  flows,  and  vor¬ 
tical  wake  flows,  the  basic  aerodynamic  characteristics  of  the  particular  problem 
are  described,  and  one  or  two  examples  of  significant  research  in  progress  are  cited. 
A  number  of  other  important  categories  and  research  issues  are  mentioned  in  passing 
in  the  final  section  of  the  paper. 


*Presented  at  the  International  Symposium  on  Recent  Advances  in  Aerodynamic  and 
Aeroacoustics,  Stanford  University,  Aug.  1983 

**Senior  Staff  Scientist,  NASA  Thermo-  and  Gas-Dynamics  Division  and  U.S.  Army 
Aeromechanics  Laboratory  (AVRADCOM) 


II.  ROTOR  BLADE  BOUNDARY  LAYERS  AND  STALL 


The  phenomenon  of  wing  stall  is  well-known  for  fixed-wing  aircraft.  It  limits 
the  minimum  flying  speed  for  takeoff  and  landing,  and  hence  determines  the  minimum 
runway  length  that  is  required  for  aircraft  operations.  On  the  other  hand,  rotor- 
blade  stall  limits  the  maximum  flight  speed  of  helicopters.  In  blade-fixed  coordi¬ 
nates,  the  difference  in  instantaneous  local  dynamic  pressure  on  the  retreating  and 
advancing  sides  of  the  rotor  disc  increases  monotonically  with  flight  speed,  result¬ 
ing  in  asymmetric  distributions  of  the  sort  illustrated  in  figure  2,  which  is  based 
on  the  experiment  of  reference  1.  To  maintain  some  semblance  of  rotor  equilibrium, 
the  blade-element  angles  of  attack  must  be  increased  on  the  retreating  blade  and 
decreased  on  the  advancing  blade.  However,  there  is  a  limit  to  how  far  the  retreat¬ 
ing  blade  angle  of  attack  can  be  increased  before  the  boundary  layer  will  separate, 
resulting  in  a  loss  of  lift  and  in  an  increase  in  blade-element  drag  and  pitching 
moment.  The  variations  in  section  lift  and  drag  translate  into  limitations  in  rotor 
thrust  and  increases  in  power  required,  respectively;  the  blade  torsion  that  is 
caused  by  the  variations  in  pitching  moment  produces  vibratory  pitch-link  loads  and 
blade  flutter. 

Retreating-blade  stall  differs  importantly  in  several  respects  from  fixed-wing 
stall  and  for  several  different  reasons,  including  the  unsteady  effects  of  rapidly- 
varying  blade-element  velocity  and  angle  of  attack,  and  the  three-dimensional  effects 
of  yaw  and  rotational  accelerations.  Experiments  indicate  that  all  of  these  factors 
may  be  important,  but  the  one  that  seems  to  be  predominant  is  the  time-dependent 
angle  of  attack,  or  pitch  rate  effect.  The  importance  of  the  rate  of  change  of  a(t) 
is  illustrated  in  figure  3,  which  shows  the  lift  and  pitching-moment  coefficient 
behavior  for  several  experimental  configurations.  All  of  the  unsteady  cases  are 
characterized  by  maximum  airloads  that  exceed  the  static  values  and  by  large  hystere¬ 
sis  in  the  flow  behavior,  according  to  whether  the  angle  of  attack  is  increasing  or 
decreasing.  The  phenomenon  of  "dynamic  stall"  is  characterized  by  strong  vortical 
disturbances  that  are  shed  from  the  leading-edge  region,  convect  over  the  upper  sur¬ 
face  of  the  airfoil,  and  induce  highly  nonlinear  fluctuating  airloads  (ref.  2).  Max¬ 
imum  lift  is  also  known  to  be  augmented  by  three-dimensional  yaw  effects  (ref.  3), 
but  the  hysteresis  evident  in  the  model  rotor  measurements  is  only  produced  by 
unsteady  effects. 

To  date,  most  of  the  research  related  to  retreating-blade  stall  has  been  per¬ 
formed  for  the  simpler  model  problem  of  a  two-dimensional  oscillating  airfoil 
(fig.  3(a)),  and  for  the  associated  unsteady  boundary  layer  behavior.  Furthermore, 
most  of  what  is  known  today  about  the  characteristics  and  various  regimes  of  dynamic 
stall  (figs.  4  and  5),  has  come  from  experiments,  in  two  or  three  dimensions. 

Attempts  to  calculate  the  quantitative  effects  of  dynamic  stall  have  not  been  very 
successful  up  to  now.  Of  the  various  two-dimensional  semi-empirical  correlation 
techniques  that  are  available  (refs.  2,  5-8),  only  the  latest  (ref.  7)  has  begun 
to  demonstrate  the  capability  of  reproducing  experimental  results  that  are  signifi¬ 
cantly  different  from  the  data  sets  that  were  used  originally  to  define  or  tune  the 
methods.  The  three-dimensional  aspects  of  the  problem  remain  virtually  untouched. 

Three  principal  issues  have  thus  emerged  in  the  general  area  of  retreating-blade 
stall.  The  first  question  is,  how  can  the  unsteady  airloads  be  predicted  with  confi¬ 
dence  over  a  wide  range  of  unsteady  flow  conditions  and  blade  geometries?  An  ade¬ 
quate  data  base  already  exists  for  the  unsteady  airloads  on  oscillating  airfoils, 
although  additional  detailed  measurements  of  boundary-layer  separation  characteris¬ 
tics  would  be  useful.  More  importantly,  a  wide  range  of  theoretical  and  numerical 
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studies  of  unsteady  viscous  flows  is  sorely  needed  to  provide  a  better  basic  under¬ 
standing  of  the  phenomena,  to  improve  the  predictive  capability,  and  to  assess  and 
guide  the  development  of  more  general  empirical  correlation  techniques. 


The  second  question  is,  to  what  extent  can  information  derived  from  either  theo¬ 
retical  or  experimental  studies  of  dynamic  stall  in  two  dimensions  be  applied  to  the 
real  three-dimensional  problem  of  retreating-blade  stall  on  helicopter  rotors?  This 
issue  can  probably  only  be  addressed  with  new  special-purpose  experiments,  on  both 
rotating  and  nonrotating  three-dimensional  models. 

Finally,  perhaps  the  most  important  issue  is  whether  practical  means  can  be 
developed  to  delay  retreating-blade  stall  or  alleviate  its  adverse  consequences. 

Since  experiments  have  shown  even  the  qualitative  features  of  dynamic  stall  to  be 
surprisingly  insensitive  to  airfoil  geometry,  major  improvements  in  stall  delay  or 
control  would  seem  to  require  a  combination  of  creative  new  ideas  for  boundary-layer 
control  or  other  unconventional  devices,  coupled  with  a  series  of  careful  experiments. 

Figures  6-8  show  three  possibilities  for  affecting  large  changes  in  dynamic 
stall  characteristics,  by  pneumatic  or  mechanical  means.  Although  these  sorts  of 
devices  have  been  used  successfully  in  fixed-wing  applications,  only  circulation  con¬ 
trol  airfoils  have  been  incorporated  into  an  experimental  helicopter  rotor  (ref.  9). 
However,  much  remains  to  be  learned  about  the  performance  of  either  leading-edge  or 
trailing-edge  blowing  under  unsteady,  three-dimensional,  and  separated  flow  condi¬ 
tions.  The  situation  is  likely  to  be  further  complicated  by  the  development  of  tran¬ 
sonic  flow  near  the  leading  edge  (see  Section  III),  if  the  blowing  delays  the  onset 
of  dynamic  stall  to  significantly  higher  values  of  lift  coefficient;  that  is,  if  it 
is  successful. 

With  regard  to  mechanical  devices,  Carr  and  McAlister  (ref.  11)  have  found  the 
leading-edge  slat  configuration  shown  in  figure  8  to  be  effective  in  suppressing 
dynamic  stall  on  an  airfoil  undergoing  sinusoidal  oscillations  in  pitch,  as  shown  in 
figure  9.  Their  particular  test  configuration  of  a  fixed -geometry  slat  suffered  from 
high  values  of  CD  at  low  angles  of  attack,  but  the  dramatic  improvement  in  Cl 
and  C^  indicates  that  the  concept  is  worthy  of  further  exploration  and  refinement. 

As  in  the  blowing  cases,  three-dimensional  and  transonic  effects  will  also  have  to  be 
studied  before  we  can  expect  this  device  to  be  incorporated  into  an  actual  rotor.  In 
any  case,  this  or  any  other  boundary-layer  control  device  that  could  enable  a  rotor 
blade  to  approach  the  high-lift  capabilities  of  modern  airplane  wings  would  pay  hand¬ 
some  dividends  for  helicopter  rotors. 


III.  TRANSONIC  AERODYNAMICS 


Transonic  flow  phenomena  commonly  appear  on  the  advancing  blade  tips  of  most 
modern  helicopter  rotors,  where  the  blade-element  Mach  numbers  are  high;  for  example, 
as  shown  in  figure  2(a).  It  should  be  emphasized  that  shock  waves  may  also  appear  on 
the  retreating  blade  as  well,  when  the  blade-element  angle  of  attack  is  sufficiently 
high  to  produce  local  regions  of  supersonic  flow  near  the  leading  edge.  Thus  the 
transonic  airfoil  problem  for  rotors  includes  both  of  the  facets  illustrated  in  fig¬ 
ure  10.  The  periodic  development  and  decay  of  shock  waves  on  rotor  blades  in  forward 
flight  lead  to  vibratory  airloads,  degradations  in  performance  and  aerodynamic  effi¬ 
ciency,  and  excessive  noise. 


In  addition  to  the  unsteady  complications,  the  transonic  flow  on  the  advancing 
blade  is  highly  three-dimensional,  generally  rendering  two-dimensional  information 
much  less  valuable  than  for  the  retreating-blade  stall  discussed  above.  Figure  11 
shows  representative  spanwise  variations  in  the  pressure  coefficient  at  the  quarter- 
chord  position  on  the  advancing  blade,  along  with  an  indication  of  the  various  flow 
regimes  that  develop.  The  change  in  the  chordwise  pressure  distribution  along  the 
span  is  also  indicated.  In  the  figure,  the  transonic  similarity  parameter,  following 
Isom  (ref.  12),  is  the  following: 


1  -  M^(r+p  sin 

K  =  -  (1) 

[(y+l)M^(r+p  sin  ip)2 i]2,i 

where  Mj  is  the  tip  Mach  number  in  hover,  =  ftt  is  the  azimuthal  position  of  the 
blade,  r  is  the  nondimens ional  radial  dimension,  p  =  V^/f.R  is  the  advance  ratio, 
t  is  the  airfoil  thickness  ratio,  and  y  is  the  ratio  of  specific  heats.  The 
inverse  of  K  essentially  represents  the  order  of  magnitude  of  the  most  important 
nonlinear  term  in  the  governing  flow  equation  relative  to  the  size  of  the  linear 
terms.  In  practice,  the  relative  importance  of  the  various  regions,  especially 
Region  2  vs  Region  4  vs  Region  5,  depends  upon  the  flight  conditions  and  the  blade 
geometry.  In  all  cases,  however,  the  flow  in  the  outboard  regions  is  three- 
dimensional,  and  calculations  also  show  important  unsteady  effects  to  occur  in  for¬ 
ward  flight,  when  the  flow  there  is  transonic. 

It  is  important  to  recognize  three  essential  differences  between  the  transonic 
aerodynamics  of  fixed-wing  aircraft  and  helicopter  blades.  The  first  is  that  the 
blade-element  Mach  number  increases  linearly  with  distance  from  the  axis  of  rotation. 
This  enhances  transonic  effects  near  the  blade  tip,  in  opposition  to  and  in  competi¬ 
tion  with  the  three-dimensional  tip-relief  effect  there.  This  tip  relief,  of  course, 
is  basically  the  same  as  for  a  nonrotating  wing.  However,  the  behavior  beyond  the 
tip  is  different,  because  of  the  linearly-varying  Mach  number;  see  reference  13. 

The  second  important  difference  is  that  in  forward  flight  the  blade-element  Mach 
number  and  angle  of  attack  vary  with  time,  or  rotor  blade  azimuthal  position,  as 
shown  in  figure  2.  This  effect  is  relatively  unimportant  in  the  inboard,  linear 
regimes,  but  the  outboard  transonic  flow  field  takes  longer  to  develop,  and  unsteady 
effects  become  essential  features  of  the  flow  at  high  advance  ratio.  Because  of  the 
lag  in  adjusting  to  the  time-  and  space-varying  local  Mach  numbers,  the  flow  in  the 
second  quadrant,  90°  <  <p  <  180°,  generally  has  stronger  and  more  persistent  shock 
waves  than  that  of  the  first  quadrant,  0  <  <p  <  90°. 

The  third  difference,  also  of  great  importance,  is  the  presence  of  concentrated 
tip  vortices  that  trail  from  the  tips  of  preceding  blades.  Depending  on  the  flight 
conditions,  these  trailing  vortices  may  remain  near  the  path  of  the  advancing  blade, 
causing  large  disturbances  to  the  shock-wave  development  and  decay.  This  aspect  of 
the  problem  cannot  be  uncoupled  from  the  rotor  wake  structure,  described  in  Section 
IV,  making  it  probably  the  most  difficult  and  challenging  aspect  of  the  aerodynamics 
of  the  advancing  blade. 

Still  another  factor  that  is  becoming  increasingly  important  is  the  trend  toward 
exotic  tip  geometries  for  advanced  rotor  designs.  Figure  12,  adapted  from  refer¬ 
ence  14,  illustrates  some  examples  which  are  already  used  on  modern  helicopters  or 
are  being  developed.  These  various  planforms  are  designed  to  alleviate  shock-wave 
development  by  sweep  effects,  to  reduce  vibratory  loads  by  either  aerodynamic  or 
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structural  dynamic  changes,  to  increase  hover  performance  by  changes  in  spanwise  cir¬ 
culation  distribution,  or  to  combine  several  of  these  factors.  The  systematic  evalu¬ 
ation  of  such  tip  shapes  is  a  difficult  task  beyond  our  present  capabilities.  In 
addition,  each  of  the  planforms  in  figure  12  would  require  some  adaptation  of  the 
computational  grid  in  any  numerical  analysis  of  the  tip  region.  Thus,  the  subject  of 
grid  generation  has  become  yet  another  fertile  area  for  research. 

The  accurate  treatment  of  the  transonic  regimes  on  rotor  blades  is,  therefore,  a 
formidable  challenge.  Fortunately,  considerable  assistance  can  be  obtained  from  the 
large  background  of  computational  aerodynamics  methods  that  have  been  developed  over 
the  past  decade  for  fixed-wing  applications,  despite  the  aforementioned  differences. 
Although  it  is  by  no  means  a  trivial  exercise  to  include  both  unsteady  and  three- 
dimensional  effects,  existing  solution  algorithms  for  the  transonic  small-disturbance, 
full  potential,  Euler,  and  even  the  Navier-Stokes  equations,  with  sufficient  effort, 
could  be  adapted  to  helicopter  problems.  Significant  progress  has  already  been  made 
in  recent  years;  and  the  level  of  complexity  that  has  been  implemented  includes  the 
following: 

(1)  The  unsteady,  three-dimensional  transonic  small-disturbance  equation  for  a 
lifting  blade  with  an  approximate  wake  model,  e.g.,  figure  13  from  refer¬ 
ence  15. 

(2)  The  quasi-steady  full-potential  equation  for  a  highly-swept,  nonlifting 
blade,  e.g.,  figure  14  from  reference  16. 

(3)  The  unsteady  full-potential  equation  for  a  nonlifting  blade,  e.g.,  figure  15 
from  unpublished  work  by  I.  C.  Chang  at  NASA  Ames  Research  Center. 

It  should  be  mentioned  that  for  realistic  rotor  calculations,  the  transonic 
fisld  calculations  have  to  be  coupled  with  the  rotor  wake,  which  is  discussed  in  more 
detail  in  Section  IV.  At  present  it  is  not  practical  to  do  more  than  solve  the  gov¬ 
erning  equations  by  finite-difference  methods  on  an  isolated  blade  and  within  a 
limited  computational  volume  around  the  blade  tip.  Therefore,  the  information  about 
the  vortex  system  in  the  wake  and  the  associated  rotor-induced  inflow  is  accounted 
for  either  as  prescribed  modifications  to  the  boundary  conditions  on  the  outer  bound¬ 
ary  of  the  computational  box,  or  as  an  additional  interior-flow  boundary  condition 
that  is  prescribed  within  the  computational  box,  or  both.  This  is  indicated  schemat¬ 
ically  in  figures  16  and  17.  As  shown  in  figure  17,  the  interior  of  the  computa¬ 
tional  domain  includes  the  usual  trailing  sheet  of  vorticity  behind  the  blade,  which 
is  assumed  to  be  undistorted,  and  one  or  more  segments  of  line  vortices  representing 
the  trailing  tip  vortex  (or  vortices)  from  the  preceding  blade(s).  The  line  vortex 
segments  are  connected  with  the  outer  boundary  by  branch-cut  surfaces,  across  which 
jumps  in  potential  are  imposed  that  correspond  to  the  strengths  of  the  concentrated 
vortices. 

The  foregoing  procedure  makes  it  possible  to  perform  meaningful  transonic  calcu¬ 
lations  for  special  cases.  However,  further  improvements  are  required  in  two  rather 
diverse  directions.  First,  the  existing  codes  are  still  too  expensive  and  complex, 
in  general,  for  routine  use  in  engineering  design.  Therefore,  serious  efforts  are 
needed  to  streamline  the  methodology  and  to  improve  the  computational  efficiency, 
without  losing  the  essential  features  of  the  problem.  On  the  other  hand,  the  current 
assumptions  and  restrictions  are  still  too  severe,  in  general,  to  allow  for  accurate 
simulations  of  many  of  the  most  crucial  operational  conditions  of  advanced  rotor 
systems.  Therefore,  improved  transonic  aerodynamics  modeling  is  needed,  but  it  must 
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not  add  excessive  burdens  to  the  computational  requirements.  Clearly  this  area 
offers  considerable  opportunities  and  challenges  for  future  research. 


IV.  ROTOR  WAKES  AND  VORTICES 


A.  Global  Wake  Modeling 

Accurate  prediction  of  the  vortical  wake  of  a  helicopter  rotor  is  probably  the 
most  important,  the  most  studied,  and  the  most  difficult  aspect  of  helicopter  aerody¬ 
namics;  consequently,  it  can  be  argued  that  this  topic  provides  one  of  the  greatest 
opportunities  ft,r  new  improvements.  Certainly  the  treatment  of  the  wake  is  crucial 
in  terms  of  performance,  efficiency,  structural  vibrations,  and  aerodynamic  noise. 
Current  methods  of  analysis  range  from  relatively  simple  momentum-theory  applications 
to  free-wake  lifting-surface  codes  that  consist  of  several  thousand  lines  of  Fortran 
statements  and  that  require  approximately  one  million  words  of  computer  memory.  In 
between  these  two  extremes,  there  are  a  variety  of  so-called  prescribed-wake  models, 
which  rely  on  some  degree  of  empiricism  to  determine  the  position  of  the  wake  vor¬ 
tices;  then  the  Biot-Savart  law  is  used  to  calculate  the  velocity  field  that  is 
induced  by  these  vortices  at  the  plane  of  the  rotor  blades. 

The  complexity  of  the  problem  in  comparison  with  fixed-wing  aircraft  is  illus¬ 
trated  in  figures  18  and  19.  The  detailed  structure  of  an  airplane  wake  is  not  all 
that  simple,  but  because  the  wake  elements  trail  rearward  along  approximately  straight 
lines  that  are  parallel  to  the  flight  direction,  simple  approximations  to  the  tip- 
vortex  rollup  and  downstream  convection  normally  suffice.  This  is  not  at  all  the  case 
for  the  helicopter  rotor,  whose  blade  tips  trace  out  prolate  cycloidal  paths  in 
space.  This  provides  numerous  opportunities  for  complex  interactions  between  the  vor¬ 
tices  and  the  blades  and  between  the  vortices  themselves. 

It  is  beyond  the  scope  and  intent  of  this  paper  to  review  the  state  of  the  art 
in  rotor  wake  modeling;  a  few  examples  are  given  to  point  out  some  features  of  the 
problem  and  to  mention  some  of  the  limitations  of  existing  methods.  For  example, 
figure  19  from  reference  17  is  representative  of  the  latest  generation  of  prescribed 
wake  analyses;  this  is  an  extension  to  forward  flight  of  one  of  the  better-known  gen¬ 
eralized  prescribed-wake  methods  for  hover  (ref.  18).  This  method  is  simpler  and 
less  expensive  than  the  free-wake  codes  of  comparable  accuracy. 

However,  an  inherent  limitation  of  the  prescribed-wake  approach  is  that  the 
empirical  determination  of  the  wake  shape  ignores  some  of  the  details  of  the  particu¬ 
lar  case  urader  consideration,  including  the  mutual  interaction  between  vortex  ele¬ 
ments.  As  a  result,  a  prescribed-wake  configuration  is  not,  in  general,  a  valid 
solution  to  an  inviscid  free-vortex  flow  (ref.  19).  In  most  cases  where  the  wake 
geometry  is  not  too  different  from  the  configurations  that  produced  the  original 
empiricism,  this  may  not  be  a  serious  error.  However,  prescribed-wake  models  are 
unreliable  for  unusual  blade  planforms  and/or  twist  distributions,  and  these  are 
often  the  cases  that  are  the  most  interesting  to  explore. 

In  some  of  these  more  challenging  cases,  the  blade  airloads  depend  strongly  upon 
both  the  wake  geometry  and  the  method  for  calculating  the  flow.  A  relatively  simple 
example  of  this  is  given  in  figure  20.  This  figure,  prepared  by  Dr.  Chee  Tung  of  the 
U.S.  Army  Aeromechanics  Laboratory,  shows  a  comparison  of  lifting-line  and  lifting- 
surface  calculations  with  the  hover  data  of  Caradonna  and  Tung  (ref.  20)  for  a  model 
rotor  with  two  low-aspect-ratio  blades  of  rectangular  planform  and  zero  twist.  The 
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computed  results  were  obtained  using  the  computer  code  described  in  reference  21. 

Here  the  local  lift  coefficient  is  nondimensionalized  by  the  tip  speed,  DR. 

For  this  example,  the  lifting-line  calculations  were  in  error,  regardless  of  the 
wake  model.  The  prescribed-wake  option  of  the  program  also  gave  an  erroneous  span- 
wise  load  distribution,  regardless  of  which  representation  was  used  for  the  surface 
of  the  blade,  primarily  because  the  prescribed  wake  geometry  was  not  correct.  On  the 
other  hand,  the  free-wake  calculation  predicted  the  actual  wake  geometry  reasonably 
well.  Consequently,  the  lifting-surface  calculation  gave  essentially  the  same  satis¬ 
factory  results  either  when  using  the  wake  prescribed  according  to  the  experimental 
measurements  or  when  using  the  free-wake  option  of  the  code.  As  a  final  footnote, 
however.  Dr.  Tung  reports  that  this  same  free-wake,  lifting-surface  code  fails  to 
predict  adequately  some  cases  with  highly  nonlinear  twist  distributions.  This  exer¬ 
cise  illustrates  the  difficulty  and  importance  of  the  wake  geometry  for  the  problem 
of  a  rotor  in  hover,  and  it  indicates  that  further  work  is  needed  even  without  the 
complications  of  forward  flight. 

B.  Blade-Vortex  Interaction 

The  strong  interaction  between  a  segment  of  a  rotor  blade  and  the  concentrated 
tip  vortices  in  the  wake  is  an  important  potential  source  of  noise  and  vibration  at 
low  and  moderate  flight  speed.  The  generic  problem,  sketched  in  figure  21,  is  an 
unsteady,  three-dimensional  close  encounter  of  a  curved  line  vortex,  at  an  arbitrary 
intersection  angle,  A^,  with  a  high-aspect-ratio  lifting  surface  that  is  executing 
combined  rotational  and  translational  motion  at  transonic  speeds.  For  many  practical 
applications,  it  is  almost  impossible  to  separate  the  transonic  aspects  of  the  prob¬ 
lem  from  the  details  of  the  vortex  structure  in  the  wake. 

Insofar  as  basic  research  opportunities  are  concerned,  the  limiting  cases  of 
A^  =  0  and  90°,  as  shown  in  figure  22,  are  of  particular  interest.  The  former  case 
is  fundamentally  unsteady  but  approximately  two-dimensional,  whereas  the  latter  is 
essentially  steady  but  highly  three-dimensional.  These  two  limiting  cases  are  the 
basis  of  several  ongoing  research  programs  that  are  designed  to  explore  the  basic 
features  of  blade-vortex  interactions,  to  develop  methods  of  calculating  such  inter¬ 
actions,  and  to  determine  the  minimum  level  of  complexity  that  will  be  required  in 
the  future  for  adequate  predictive  capability  in  the  complete  rotor  environment.  The 
case  of  a  vortex  that  is  parallel  to  the  free  stream  has  received  more  attention  in 
the  past,  as  it  is  more  relevant  to  fixed-wing  applications  and  it  is  relatively  easy 
to  set  up  experimentally.  On  the  other  hand,  the  case  of  the  vortex  that  is  parallel 
to  the  leading  edge  is  a  simpler  problem  for  theoretical  and  numerical  analysis,  and 
it  requires  less  computer  memory  to  store  the  solution  variables. 

For  illustrative  purposes,  we  shall  briefly  review  some  recent  numerical  work 
for  the  two-dimensional  unsteady  case,  A^  ■  0.  Two  of  the  basic  issues  are  the  intro¬ 
duction  of  the  vortex  into  the  numerical  computations  and  the  determination  of  its 
effect  on  the  unsteady  pressure  distribution  on  the  airfoil  and  nearby.  The  vortex 
can  be  introduced  in  two  different  ways,  as  indicated  in  figure  23.  Caradonna,  et  al 
(ref.  15)  and  George  and  Chang  (ref.  22)  have  utilized  the  branch-cut  method,  as 
described  in  Section  III,  and  McCroskey  and  Goorjian  (ref.  23)  employed  the 
prescribed-disturbance  method.  All  three  groups  solved  the  transonic  small- 
disturbance  equations  by  time-accurate  ADI  methods  that  had  been  developed  and 
checked  out  previously,  and  the  various  results  are  in  reasonable  agreement  in  most 
cases. 
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Representative  results  from  reference  23  are  shown  in  figures  24  -  25  for  two 
symmetrical  airfoil  sections.  Figure  24  shows  the  distortions  in  the  chordwise  pres¬ 
sure  distributions  on  the  airfoil  as  the  vortex  passes  underneath.  Similar  results 
are  shown  in  figure  25  for  a  thinner  airfoil,  including  the  time  histories  at  indivi¬ 
dual  locations  on  the  airfoil.  These  results  illustrate  that  strong  gradients  in 
pressure  occur  with  respect  to  both  time  and  space,  due  to  vortex  encounter.  These 
gradients  can  be  especially  significant  in  the  leading-edge  region  of  a  thin  airfoil. 

Although  the  subject  of  airfoil-vortex  interaction  is  just  beginning  to  receive 
significant  attention,  several  important  features  of  the  problem  have  already  emerged. 
First,  the  effect  of  the  vortex  is  felt  primarily  through  the  vertical  velocity  that 
it  induces,  which  to  first  order  appears  as  a  time-dependent  perturbation  in  the 
effective  angle  of  attack,  and,  secondarily,  through  its  horizontal  induced  velocity. 
Second,  unsteady  lag  effects  are  very  important,  especially  in  the  transonic  case; 
calculations  show  enormous  differences  between  quasi-steady  and  unsteady  solutions 
and  between  the  results  for  the  vortex  locations  upstream  and  downstream  of  the  air¬ 
foil.  Third,  other  strong  nonlinear  effects  have  been  noted;  even  the  qualitative 
pressure  variations  with  respect  to  time  and  space  differ  markedly  from  the  predic¬ 
tions  of  linear  theory,  and  the  vortex  distorts  the  flow  on  the  nearest  surface  of 
the  body  much  more  than  on  the  opposite  surface.  Finally,  some  of  the  pressure  per¬ 
turbations  appear  to  leave  the  body  as  radiating  waves;  this  phenomenon  would  seem  to 
represent  the  far-field  blade-vortex  interaction  noise  of  helicopter  rotors. 

The  studies  to  date  (refs.  15,  22-23)  have  been  performed  with  the  small- 
disturbance  equations,  whose  accuracy  is  questionable  in  the  very  region  that  seems 
to  be  the  most  important  (ref.  13);  that  is,  near  the  leading  edge.  Therefore,  it 
seems  imperative  to  check  the  small-disturbance  codes  with  more  accurate  formulations 
and  with  experiments.  Also,  an  ideal  potential  vortex  with  an  invariant  structure 
has  been  assumed  thus  far;  but  in  reality,  close  encounters  probably  alter  the  vortex 
core  significantly  and  may  lead  to  vortex  bursting.  For  rotor  applications,  future 
studies  must  be  extended  to  include  three-dimensions,  where  it  will  be  even  more 
important  to  establish  what  minimum  level  of  complexity  in  the  governing  equations 
will  suffice  and  to  determine  the  most  expeditious  way  to  introducing  the  vortex  into 
the  computational  domain.  Finally,  the  severity  of  the  blade-vortex  interaction 
appears  to  be  highly  sensitive  to  the  strength  and  position  of  the  vortex,  especiaxjLy 
for  transonic  flow;  therefore,  the  results  of  the  computations  will  be  no  better  than 
the  predictions  of  the  wake  model  used. 


V.  CONCLUDING  REMARKS 


The  examples  discussed  in  the  preceding  sections  are  but  a  few  of  the  myriad  of 
serious  aerodynamic  problems  of  helicopters  that  provide  rich  opportunities  for 
future  research.  There  are  other  examples  which  are  not  discussed  in  this  selective 
review;  however,  there  are  several  that  are  worth  mentioning  in  passing. 

The  first  is  the  general  category  of  aerodynamic  interference.  As  indicated  in 
figure  1,  the  flow  field  of  the  various  components  of  the  helicopter  frequently 
influence  other  flow  regions,  and  these  interactions  can  be  complex  and  nonlinear. 
Typical  examples  include  the  mutual  interference  and  interaction  between  the  main 
rotor  and  the  fuselage,  between  the  main  rotor  and  the  tail  rotor,  between  the  rotor 
hub  and  the  fuselage,  and  between  the  tail  rotor  and  the  tail  surfaces.  Panel 
methods  are  being  developed  for  predicting  these  flows;  e.g.,  reference  24,  but  the 
curient  industry  approach  is  mostly  empirical.  A  closely-related  problem  that  has 
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arisen  recently  is  the  rotor-induced  ’’download"  on  the  wing  of  the  tilt-rotor  air¬ 
craft,  figure  26.  In  hover,  the  downwash  from  the  rotors  impinges  on  the  wing,  caus¬ 
ing  a  net  loss  in  lifting  capability  of  10%  or  more  of  the  gross  weight  of  the  air¬ 
craft.  Predicting  this  effect  accurately  and  designing  flap  configurations  to 
minimize  the  download  are  particularly  challenging  problems,  but  ones  which  have 
great  payoff. 

Another  important  area  is  that  of  simply  improving  the  aerodynamic  tools  that 
are  available  to  helicopter  engineers.  The  need  for  simpler  and  more  efficient  compu¬ 
tational  techniques  that  can  more  accurately  solve  more  complex  problems  has  already 
been  stressed.  However,  improvements  in  the  coupling  of  aerodynamic  prediction  tech¬ 
niques  with  structural-dynamics  codes  is  also  essential.  In  addition,  there  would 
appear  to  be  tremendous  opportunities  for  coupling  efficient  aerodynamic  codes  with 
numerical  optimization  techniques  to  design  improved  blade  configurations.  This  has 
already  been  demonstrated  in  two  dimensions  for  transonic  rotor  airfoil  sections 
(ref.  25).  The  extension  to  three  dimensions,  that  is,  including  planform  and  twist 
distributions,  is  the  next  logical  step,  although  it  is  obviously  a  very  large  step. 

Experiments  will  continue  to  play  a  major  role  in  helicopter  technology  for  the 
foreseeable  future,  and  there  is  much  to  be  done  in  the  area  of  instrumentation, 
experimental  methods,  and  testing  techniques.  One  perennial  issue  is  the  gap  between 
model  and  full-scale  testing;  the  recent  paper  by  Keys,  et  al  (ref.  26),  for  example, 
discusses  some  aspects  of  this  problem.  Mach  number  scaling  is  clearly  essential; 
but  this  is  hardly  a  research  issue  any  more,  even  though  constructing  dynamically- 
scaled,  instrumented  model  rotors  to  operate  at  full-scale  tip  speeds  can  be  a  vexing 
engineering  exercise.  On  the  other  hand,  typical  Reynolds  numbers  for  model  and  full- 
scale  helicopters  span  regimes  where  the  boundary-layer  transition  and  separation 
characteristics  may  be  considerably  different.  Therefore,  the  model/full-scale  issue 
will  continue  to  be  important  for  any  aspect  of  helicopter  aerodynamics  in  which 
viscous  effects  play  a  significant  role. 

A  strong  need  for  a  wide  variety  of  special-purpose  experiments  will  continue  to 
exist.  For  example,  tests  of  complete  rotor  systems  certainly  have  their  place,  but 
a  full-scale  helicopter  blade  is  not  the  best  experimental  configuration  to  study  the 
effects  of  unsteadiness  on  the  Reynolds  stresses  of  a  separating  turbulent  boundary 
layer,  nor  to  validate  initially  a  new  transonic  code.  Relatively  simple  experiments 
that  are  run  under  well-controlled  and  thoroughly-documented  conditions  are  essential 
to  sort  out  the  dominant  features  of  the  various  flow  fields,  to  increase  our  funda¬ 
mental  understanding  of  the  basic  phenomena,  to  guide  the  development  of  new  theoret¬ 
ical  models,  and  to  validate  methodically  the  new  computational  algorithms  and 
prediction  techniques  as  they  are  being  developed. 

Finally,  for  both  the  global  studies  and  for  the  special-purpose  experiments,  a 
wide  variety  of  new  instrumentation,  measurement  techniques,  and  data  acquisition  and 
processing  methods  need  to  be  developed.  These  are  required  to  cope  with  the  com¬ 
plexity  of  the  unsteady  three-dimensional  flow  fields  of  the  various  helicopter  com¬ 
ponents  and  with  handling  the  large  volume  of  data  that  are  generated.  It  is  as 
challenging  to  measure  the  right  information  and  to  eliminate  the  superfluous,  as  it 
is  to  develop  the  new  predictive  methodology. 

In  conclusion,  the  examples  discussed  in  this  paper  illustrate  the  complexity  of 
helicopter  aerodynamics  and  the  many  important  issues  that  remain  to  be  studied. 

There  are  numerous  challenges  and  a  wide  range  of  opportunities  for  the  technologist 
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and  for  the  basic  research  scientist  to  explore.  Their  various  contributions  will 
assure  continuing  major  improvements  to  this  unique  flying  machine. 
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(b)  Unsteady  yawed  flow,  pitching  oscillations  (ref.  3) 


Figure  3.-  Continued. 
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(c)  Model  rotor  blade  in  forward  flight  (ref.  4) 
Figure  3.-  Concluded. 


Figure  4.-  Light  dynamic  stall  on  an  oscillating  airfoil 
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Figure  5.-  Deep  dynamic  stall  on  an  oscillating  airfoil 
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Figure  7.-  Leading-edge  blowing  to  alleviate  retreating  blade  stall  (ref.  10) 


Figure  8.-  Leading-edge  slat  to  alleviate  retreating  blade  stall  (ref.  11). 
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Figure  11.-  Spanwise  variations  and  flow  regimes  for  an  advancing  rotor  blade. 
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Figure  12.-  Examples  of  advanced  rotor  blade  tips. 
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Figure  13.-  Pressure  distributions  on  a  lifting  rotor  at  high  advance  ratio,  using  an 
unsteady  small-disturbance  code  (ref.  15);  u  =  0.36,  CT/o  =  0.05,  <J>  =  90°. 
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Figure  14.-  Pressure  distributions  on  a  nonlifting  rotor  at  high  advance  ratio,  using 
a  quasi-steady  full-potential  code  (ref.  16);  u  ■  0.40,  i|»  ■  90°. 


22 


-  UNSTEADY 

- QUASI-STEADY 


Figure  15.-  Pressure  distributions  on  a  nonlifting  rotor  at  high  advance  ratio,  using 
an  unsteady  full-potential  code;  p  =  0.55,  ip  =  120°. 
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Figure  16.-  Aerodynamic  modeling  of  the  rotor  for  finite-difference  methods. 


Figure  17.-  Boundary  conditions  for  potential  flow  past  rotor  blades. 


Figure  18.-  Wake  of  a  rotor  in  hover  vs.  a  fixed  wing  wake 
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Figure  19.-  Wake  of  a  rotor  in  forward  flight  (ref.  17):  4  blades,  p  =  0.1 


Figure  20.-  Concluded 
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Figure  21.-  Helicopter  blade-vortex  interaction. 
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Figure  22.-  Limiting  cases  of  blade-vortex  interaction. 
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Figure  23.-  Two  methods  of  introducing  a  vortex  into  the  flow  field 
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Figure  24.-  Transonic  airfoil-vortex  interaction  (ref.  23) 
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(a)  Time  histories  at  specific  points  on  the  airfoil. 
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(b)  Pressure  distributions  on  the  airfoil  at  specific  times. 

Figure  25.-  Pressure  variations  on  the  lower  surface  during  an  airfoil 
interaction  (ref.  23)  NACA  64A006  airfoil,  MOT  =  0.85,  a=0,  r=0. 
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Figure  26.-  The  tilt-rotor  aircraft  in  the  hover  mode 
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